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Abstract. We briefly review the theory for electromagnetic reactions in light nuclei based 
on the coupled-cluster formulation of the Lorentz integral transform method. Results on 
photodisintegration reactions of and are reported on and preliminary calcula¬ 
tions on the Coulomb sum rule for ''He are discussed. 


1 Introduction 

The investigation of nuclear reactions from first principles is fundamental in bridging nuclear physics 
with the underlying quantum chromo-dynamics regime. Nuclear reactions induced by electromagnetic 
probes are very useful as the electromagnetic current is well known and a clean comparison with 
experimental data can be performed Hi. Nowadays this valuable information is not only accessible 
for the lightest nuclei, but novel theoretical approaches are being developed to tackle nuclei with a 
larger number of nucleons. Below, we briefly review the Lorentz integral transform (LIT) approach in 
its coupled-cluster (CC) theory formulation Q. Thereafter, we discuss some recent results regarding 
photoabsorption and electron-scattering reactions. 


2 Theoretical Formulation 

The key ingredient to study reactions induced by electromagnetic external probes, like photons or 
electrons, is the nuclear response function 

R(cj,q) - ^ \{n\ 0(q) \0)f6(E, - Eq - cu) , (1) 

n 

where Oiq) is the excitation operator, which will depend specifically on the external probe and on 
the momentum-transfer q. The nuclear response function is a dynamical observables which requires 
knowledge on the whole spectrum of the nucleus, being |0) and \n) the ground- and excited-state, 
respectively. The exact solutions of the many-body problem for the excited-states, typically in the 
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continuum, is limited in mass number and in energy u). In particular, it is difficult to calculate when 
many channels are open. However, one can use the LIT approach ||2l to reduce the continuum problem 
to the solution of a bound-state equation. In particular, recently this method was reformulated using 
coupled-cluster theory Q, obtaining a new approach to extend ab-initio studies of electromagnetic 
observables to heavier and neutron-rich nuclei EH. 

In coupled-cluster theory one introduces the similarity transformed Hamiltonian 

H = exp(-7’)Hwexp(r), (2) 

where Hff is normal-ordered with respect to a chosen Slater determinant |<l)o). Correlations are intro¬ 
duced through the cluster operator T which is expanded into particle-hole {ph) excitation operators, 
i.e.,T - T\ + T 2 + ..with the \p-\h excitation operator T\, the 2p-2h excitation operator T 2 , and so 
on. In the coupled-cluster formalism ||4| the response function in Eq. Q becomes 

R(to, q) = ^<OL|0^(?)|nfi)<nL|0(g')|OR)d(£'„ - Eq - w), (3) 

n 

where = e~^0^^\q)e^ is the similarity transformed (adjoint) excitation operator and (0 l|, 

|0«) ((hlI, Ihr)) are the left and right reference ground-states (excited-states), respectively (see also 
Ref. 161). The utilization of the LIT method requires to find the solution '!'« of the following bound- 
state equation EH 

(H-A£o)I'1'r) = 0(?)|O«), (4) 

where ALq is the ground-state correlation energy and |0s) = |*l>o)- By assuming that |'Tfi) = ^|<5o) 
and expanding R. in ph excitation operators, similarly to what done with T, one can see that Eq. 0 
is equivalent to an equation-of-motion E with a source in the right-hand-side. In Ref. EIH we have 
implemented this method with a truncation of both the operators T and R. up to the 2p-2h excitation 
level, namely up to single and double (SD) excitations. Such a scheme is being named here the 
LIT-CCSD method. 

By solving Eq. 0 one finds the LIT of the response function, integrated from threshold to 
infinity 

r r“ R{u),q) - - 

L(wo,r) = - day- - 2 r 2 (5) 

where cuq and T are auxiliary parameters and (T'^.l is the solution of the equivalent to Eq. 0 for 
the left-hand-side. To obtain the response function R{aj, q) we then invert Eq. Q using the method 
outlined in Refs. Eia, which looks for the regularized solution of the integral transform equation. 

Below we highlight some of our recent calculations obtained for photoabsorption and electron 
scattering reactions using as only ingredient an Hamiltonian H that contains, besides kinetic energy, 
a two-nucleon force derived in chiral effective field theory at next-to-next-to-next-to-next-to-leading 
(N^LO) order HI]. 


3 Photodisintegration reactions 

Photodisintegration reactions have been widely studied in the ’70s with experiments on a variety of 
stable nuclei. The main observed feature of the measured cross sections was a very pronounced peak, 
referred to as the giant dipole resonance, located at excitation energies of about 10-30 MeV. Eirst the¬ 
oretical interpretations were given in terms of collective models Em and until very recently, most 
microscopic calculations were obtained in terms of phenomenological approaches and interactions. 
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Figure 1. photodisintegration cross section 

compared with data from Ref. El. The curve is 
calculated with the N^LO nucleon-nucleon force 
and is shifted to the experimental threshold. 



(0[MeV] 

Figure 2. "^°Ca photodisintegration cross section 
compared to data from Ref. El- The curve is cal¬ 
culated with the N^LO nucleon-nucleon force and is 
shifted to the experimental threshold. 


see, e.g., Ref. IH. The LIT-CCSD approach introduced above, offers the opportunity to investigate 
photodisintegration reactions from an ab-initio point of view using nucleon-nucleon potentials, which 
reproduce two-nucleon scattering data IfTTl . Moreover, it allows one to study equivalent processes 
in neutron-rich nuclei, which have been more recently investigated with Coulomb excitations experi¬ 
ments at the rare isotope facilities, see e.g. Ref. ca. 

In the long-wavelength limit, the photodisintegration cross section can be written as 

o-y(oj) - 47r^a at R^{at ), (6) 


where a is the electromagnetic coupling constant and R^(at) is the dipole response function, namely 
Eq.Q with the the dipole operator as excitation operator 


D, 


Z 


"^cm) 




(7) 


Here, Zk is the z-coordinate of the k-th nucleon in the lab-frame, while Zc„, denotes the z-component 
of the center of mass of the nucleus and is the third component of the isospin of the k-\h nucleon. 

In Refs. |I3]|4| we implemented the one-body operator of Eq. (0 in the solution of Eq. Q and in¬ 
vestigated the photodisintegration of a variety of light and medium-mass nuclei. After having bench- 
marked the LIT-CCSD method with exact hyperspherical harmonics ini on '^He, where we have 
observed that higher-order corrections to the T and R operator are very small, we have proceeded to 
study heavier nuclei. In particularly, here we review our results for the neutron-rich isotope in 
Eig.[^and the stable ^°Ca nucleus in Eig.[^ 

Eor one notices the appearance of a small peak at low energy, also named pigmy dipole 
resonance and experimentally observed in a variety of neutron-rich nuclei, which agrees very nicely 
with data from Leistenschneider et al. ini. It is worth noticing that nothing in the interaction has 
been adjusted to and that no cluster structure has been imposed a priori in the calculation. This 
is not the first time that the LIT approach suggests the existence of a low-energy dipole mode, as in 
Ref. CD a similar situation was found in ®He. 
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Figure 3. LIT-CCSD sum rules from CJJ and CJ 
for '^He for different model space sizes with 
the N^LO nucleon-nucleon chiral interaction. 



q [MeV/c] 

Figure 4. Recursive sum of isovector CSR for 
'^He calculated using the chiral N^LO nucleon- 
nucleon potential in the LIT-CCSD scheme. 


Because the computational cost of the coupled-cluster method scales mildly with respect to the 
mass number, we can extend our studies to the photonuclear cross-section of '^^Ca. Fig. shows a 
comparison of the LIT-CCSD calculations with data by Ahrens et al. m. We clearly see a giant 
resonance, even though the theoretical prediction is slightly broader and lower in strength. Both in 
Fig. r and Fig. the band in the curve is obtained by inverting LIT with different F parameters in 
Eq. (5 1 and can be seen as a lower estimate of the theoretical error bar. 


4 Electron-scattering reactions: Coulomb sum rule 

Electron scattering reactions have been already investigated within the LIT method, see, e.g., uni. 
With the objective in mind to eventually tackle neutrino-nucleus interactions relevant for long-baseline 
neutrino experiments, we initiated a study of electron-scattering reactions within coupled-cluster the¬ 
ory. The latter has in fact the potential to be applied to and other medium-mass nuclei, which 
are used as detector materials in neutrino experiments. Eollowing the work done with the Green’s 
function Monte Carlo method in Ref ifTSll . we will first investigate the Coulomb sum rule. 

The longitudinal Coulomb sum rule (CSR) is defined as 


Sdq) 


1 r 
zXi 


do) 


q) 

cihQ^) 


(8) 


where Ri{cx),q) is the longitudinal response function, i.e., Eq. 0 with 0{q) - piq) and a recoil term 
^ in the delta function, while is the proton electric form factor as a function of the square of 

the four-momentum Q^. The charge operator can be written as 


piq) = ^ 






(e") 






(9) 


with G'^(Q^) being the neutron electric form factor. Similarly to what done in Ref. IfTTlI we perform a 
multipole expansion of Eq. Q in terms of Coulomb multipoles, separating isoscalar (i) and isovector 
(v) contributions as 

J 


(10) 
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Figure 5. Preliminary results for the “^He CSR from coupled-cluster theory compared with experimental CSR 
obtained from data from Buki et. al. (red triangle) ioi and from Ref. Q9) (green square), as described in the 
text. 


Each of the Coulomb multipole is calculated separately solving Eq. Q for various fixed values of the 
momentum transfer q. Below, we show LIT-CCSD results for "^He obtained using the N^LO nucleon- 
nucleon chiral interaction. 

In Fig. we test the convergence of our calculations in terms of the model space size N„uix - 
2n - 1 - /, determined by the number of single-particle shells used in the ph expansion. For the isovector 
monopole Cg and isovector dipole C\ operator, one observes that very good convergence is reached 
already for Nmax - 14. For higher order multipoles, slightly larger model spaces are needed. With 
respect to exact hyperspherical harmonics m, Cj calculations from coupled-cluster theory agree 
very nicely up to about q - 200 MeV/c, but are 20% smaller at higher momenta. Note that at energies 
of about 300 MeV and higher, the chiral interaction of Ref. HID is at the limit of its applicability. 
Work is in progress to understand the role of center of mass contamination in our calculation, given 
that the operator in Eq. (|^ is not translational invariant. 

In Fig. 1^ we compare the recursive sum of the isovector multipoles that make up the CSR. One 
can readily see that, while at momentum transfer below ^ = 100 MeV/c two or three multipoles 
suffice, at the highest momentum transfer considered q - 500 MeV/c, nine multipoles need to be 
calculated before reaching convergence. A similar situation is found for the isoscalar multipoles. 
Finally, in Fig. we compare the total CSR calculated from coupled-cluster theory with experimental 
data. Our calculation still misses the contribution of the isoscalar monopole term. Measurements 
of the longitudinal response functions at intermediate momentum transfer have been performed in 
the past and are collected in Ref. d, while low momentum data q — 200 and 250 MeV/c are 
taken from Buki et. al. ESI. Since finite maximal values of the energy transfer cOmax measured in 
experimental data, the experimental CSR is obtained as 


Sdq) 


_ 1 r 
“ z L 


• ^ Rdo), q) ^ ^ 

ciW) 


( 11 ) 


where S l. tail is taken from the theoretical calculations of ^He response functions of Ref. El. 

As shown in Fig. our calculation agrees well with the experimental data from low- to intermediate- 
momenta. It is known from Ref. IfTTlI that three-body forces are not very pronounced on the CSR at the 
values of momentum transfer where data are available. Work is being presently directed to benchmark 
these studies with hypersherical harmonics, with the prospect of tackling '®0. 
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5 Conclusions 

After a short review of the coupled-cluster formulation of LIT method, we report on our recent calcu¬ 
lations of the photodisintegration of the neutron-rich nucleus and the stable "^°Ca medium-mass 
nucleus. We then show our preliminary calculation of the Coulomb sum rule for '^He, which describe 
rather well the available experimental data. A benchmark with data is essential for any theory that is 
going to be used to model the neutrino-nucleus interaction taking place within the detectors of neu¬ 
trino long-baseline experiments. Using the LIT-CCSD theory, our long term goal is to study neutrino 
scattering off '®0, which is relevant for the T2K experiment. 
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